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We study features in tlie optically detected magnetic resonance (ODMR) signals associated with 
negatively charged nitrogen- vacancy (NV~) centers coupled to other paramagnetic impurities in 
diamond. Our results are important for understanding ODMR line shapes and for optimization 
of devices based on NV~ centers. We address mislabeling of side features to the unperturbed 
NV~ magnetic resonance. Side resonances separated by around 130 MHz are due to hyperfine 
coupling between NV~ centers and nearest-neighbor ^'^C nuclear spins. Side resonances separated 
by approximately {40, 260, 300} MHz are found to originate from simultaneous spin flipping of NV~ 
centers and single substitutional nitrogen atoms. All results are in agreement with the presented 
theoretical calculations. 

PACS numbers: 61.72.jn, 81.05.ug, 76.70.Hb 



I. INTRODUCTION 

Negatively charged nitrogen- vacancy (NV~ ) color cen- 
ters in diamond are promising candidates for magnetom- 
etry with an unprecedented combination of sensitivity 
and spatial resolutiorP. In particular, diamonds with 
high density of NV~ centers are appealing for ensemble 
magnetometrjEl. Optically detected magnetic resonance 
(ODMR) signals enable probing the energy levels of NV~ 
centers, from which the external magnetic field can be 
inferred. In this work, we analyze lineshapes in ODMR 
signals of NV~ centers in diamond. 

In the absence of magnetic fields, the NV^ center has 
a magnetic resonance at a frequency of approximately 
2870 MHz, which corresponds to a transition between 
the triplet ground-state magnetic sublevels with electron 
spin projections nig = and mg = ±1 [see Fig. l(a)| . 
Side resonances around this central resonance have been 
reported in the literature. These extra resonances are 
due to coupling between NV~ centers and other param- 
agnetic impurities in the diamond lattice. 

A pair of side resonances asymmetrically offset from 
the central peak and separated from each other by 130 
MHz have been attributed to hyperfine interaction be- 
tween a NV~ center and a nearest-neighbor ^■^C nuclear 



spirF"^. Hyperfine coupling with '^C spins located in 
other lattice sites has also been studied in detail^. Such 
interactions have been used to demonstrate quantum in- 
formation processing using NV~ centers and "'^'^C nuclear 
spinPm. 

Other side resonances have been observed in nitrogen 
rich diamond^^ These resonances become particu- 
larly pronounced at high microwave powers, where the 
central resonance saturates, while these side resonances 
continue to grow with increasing microwave power. Un- 
derstanding this regime is important for understanding 
light-narrowing effects in ODMR and for optimizing the 



performance of sensors based on ensembles of NV^ cen- 
ters in diamond—. 

At zero magnetic field, van Oort et al^ report pairs 
of side resonances separated by 140 MHz and 280 MHz, 
symmetrically displaced from the central resonance. Van 
Oort et al]^ attribute these resonances to simultaneous 
spin fiips of NV~ centers and single substitutional nitro- 
gen atoms (PI centers). However, a part of the interpre- 
tation and theoretical analysis of Ref. [TSjis questionable, 
and the resonances with reported separation of 140 MHz 
seem to be blended with the -'^^C related resonances sep- 
arated by 130 MHz. 

In this paper we present improved ODMR spec- 
tra showing the different side resonances in diamonds 
with both high and low nitrogen concentration. The 
microwave-power and magnetic-field dependences of 
these features are studied in detail. We are able to un- 
ambiguously ascribe the origins of all the observed side 
resonances. We argue that van Oort et al^ and later 
Baranov et al^ and Babunts et al^have misinterpreted 
certain side resonances. Furthermore, we show that the 
side resonances separated by 280 MHz observed in Ref.fT^ 
at zero magnetic field are partially unresolved pairs of res- 
onances; one group separated by 260 MHz and another 
by 300 MHz. We determine the origin of side resonances 
separated by approximately {40, 260, 300} MHz to be si- 
multaneous spin flips of the NV~ centers and PI centers, 
and we clearly illustrate that the origin of the resonances 
separated by 130 MHz is due to hyperfine coupling to 
^^C. These conclusions are founded upon experimental 
evidence and theoretical calculations. 



II. METHODS 

In our experiment, we used a conventional confocal- 
microscopy set-up shown in Fig. l(b)| We focus 
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FIG. 1: (a) Energy-level structure of the NV~ center. Gray 
areas represent phonon sidebands. ■^A2 and ''E are the spin- 
triplet ground and excited states, respectively. The elec- 
tron spin-projection along the NV~ axis is denoted mg. (b) 
Confocal-microscopy set-up. 



TABLE I: Fit parameters obtained from a fit of the measured 
contrast as a function of resonant microwa ve Ra bi frequency 
to Eq. ([l]). Data and fits are shown in Fig. 2(b) 



Resonance 


K 


Qsat (MHz) 


Group A 
Group C 
Group D 
Group E 


79(14) 
3.4(6) 
0.012(6) 
0.011(6) 


0.039(7) 
0.032(3) 
2.4(1) 
2.5(1) 



that the concentration of NV~ centers is 25-100 times 
larger in the W5 sample compared to the W2 sample. 
The concentration of NV~ centers in each sample varies 
by a factor of 2-4 depending on the exact spot on the 
diamond. 



III. RESULTS AND DISCUSSION 



continuous- wave 532 nm laser light onto the diamond sur- 
face and optically pump the NV~ centers to the \ms = 0) 
ground-state sublevel. The fluorescence (wavelength « 
638-800 nm) propagates back though the lens, through 
a dichroic mirror, and onto a detector (photodiodc or 
photomultiplier tube). Microwaves are applied with a 
wire pressed across the diamond surface. The microwave 
frequency is scanned over the transition between the 
\ms — 0) and the \mg — ±1) sublevels. On resonance, 
the NV~ centers are transferred to the \ms — ±1) states, 
which leads to lower red fluorescence intensity than when 
the centers are in the \ms = 0) state (see the discussion 
in Ref. [T^. We also measured Rabi oscillations by op- 
tically pumping NV~ centers to the \ms = 0) state and 
applying resonant-microwave frequency lasting a few mi- 
croseconds while keeping the laser on. This allows us 
to observe population oscillations between the \mg = 0) 
and \ms = ±1) states in real-time by monitoring fluores- 
cence intensity modulation while microwaves are on. We 
used low laser power in order to avoid optical pumping 
during the duration of oscillation. 

We tested diamond samples with different amounts 
of impurities. Single-crystal diamonds were obtained 
from Element Six with nitrogen concentrations of 
100 ppm synthesized by high-pressure, high-tcmperature 
technique (HPHT) or with nitrogen concentrations of ^ 1 
ppm synthesized by chemical-vapor deposition (CVD)^^. 
Then, the samples were irradiated with relativistic elec- 
trons to create vacancies. In the measurements described 
here, we used the HPHT W5 sample and the CVD W2 
sample, which were subjected to irradiation doses of 10^^ 
cm~^ and 10^^ cm~^, respectively. Both irradiation pro- 
cesses were performed at the Mainz Microtron facility in 
Germany. Finally, the diamond samples were annealed 
in vacuum for 3 hours at 750° C to combine vacancies 
with nitrogen atoms to complete the formation of NV~ 
centers. Based on measurements of the amount of fluores- 
cence emitted from the the two samples, we can estimate 



Figure 2(a) shows experimentally obtained ODMR 
spectra at zero magnetic field in the frequency range 
2600 MHz to 3150 MHz. At 2870 MHz, we observe 
the central resonance, which is labeled as resonance A 
in Fig. 2(a) Resonance A corresponds to a transition 
between the \ms = 0) and \mg = ±1) states of the NV^ 
center. In measurements taken at small resonant mi- 
crowave Rabi frequencies (yiR), which correspond to low 
microwave powers, the central resonance is split into two 
resonances because the presence of strain causes mix- 
ing between the \ms — ±1) sublevels and energy split- 
ting between the resulting eigenstates [see, for exam- 
ple, the measurements taken at fij? = 0.03 MHz and 
^lu = 0.30 MHz in Fig. 2(a) . Contrast, which is de- 



fined as the relative change of fluorescence intensity when 
the microwaves are on and off resonance with a transi- 
tion, increases with increasing resonant microwave Rabi 
frequency. Figure |2(b)| shows the contrast of each reso- 
nance as a function of rin. The measured contrasts for 
each resonance were fitted to the function^'^ 



1 



(1) 



where C(r2fl,) is the contrast, K is a proportionality con- 
stant, and flsat is the saturation Rabi frequency. The 
measured contrasts and the fits are shown in Fig. |2(b)[ 
and the fitted parameters are shown in Table |T] For reso- 
nance A, the saturation Rabi frequency was found to be 
0.039(7) MHz. 

Apart from this central resonance, there are two small 
resonances in group C, Ac — 126(2) MHz, asymmet- 
rically displaced to the low- and high- frequency sides 
of the center by around 56 MHz and 70 MHz, respec- 
tively. For the lowest microwave power corresponding to 
ilu — 0-03 MHz, only side-resonance group C and the 
central resonance are observed [see Fig. [2(a)] . As seen 
in Fig. |2(b)| and from the fitted parameters in Table [ij 
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FIG. 3: Mixed energy-level diagrams for (a) hyperfine cou- 
pling of NV~ centers to ^"^C nuclear spins and (b) simulta- 
neous spin flipping of NV~ centers and PI centers. Each 
transition causing a resonance in Fig. |2(a)"] is labeled with the 
appropriate letter. 



FIG. 2: (a) Side-resonance features in ODMR are observed on 
either side of the central resonance at 2870 MHz. Measure- 
ments were done at zero magnetic field with a laser power of 
0.7 mW, on sample W5. Q,r denotes the resonant microwave 
Rabi frequency. The features in this signal are too broad to 
resolve side-resonance group B. (b) Measurement contrast of 
side resonances A, C, D, and E is plotted as a function of 
57 ii, which is proportional to the square root of microwave 
power applied to the sample. The contrast of side-resonance 
group C is multiplied by 20 for clarity. We only plot the con- 
trasts of the C resonances for Rabi frequencies f2_H < 0.6 MHz, 
because the C resonances overlap with the B resonances for 
higher resonant microwave Rabi frequencies. 



the contrasts of resonance groups A and C saturate at 
similar values of the resonant microwave Rabi frequency. 
This similar growth of contrasts of side-resonance groups 
A and C is expected assuming that the C resonances 
are related to those NV~ centers which have a nearest- 
neighbor ^^C atom. The NV~ centers (which has elec- 
tronic spin 1) and ^^C atoms (which has nuclear spin 1/2) 
can interact through the hyperfine interaction, which is 
strongest when the ^^C atom is in a nearest-neighbor 
postion to the vacancy, where the electron density is 



the highest. The relative sizes of the contrasts can be 
estimated from the ratio of NV^ centers with nearest- 
neighbor ^'^C atoms to those without nearest-neighbor 
^^C atoms. Since the natural abundance of ^^C atoms 
is 1.1%, and there are three nearest-neighbor positions 
to the vacancy, we expect the ratio of contrasts to be 
roughly (3*0.011/2) = 1.65% w 1/60. The factor of 2 ap- 
pears because the ground state \ms = 0) splits into two 
states: |m g = , mj = +1/2) and \ms = Q,mi = —1/2) 
[as in Fig. 3(a)| due to the hyperfine interaction (here 
ms is the NV~ electronic spin projection and m/ is the 
^^C nuclear spin projection). In the case that the ^"^C 
nuclear spins are unpolarized, the two ground states con- 
tain half of the population each. The measured ratio of 
the contrasts is found to be 3.0(3)%, which is a bit higher 
than the expected value because the side resonances are 
less affected by strain than the central resonance. The 
reason these side resonances are less affected by strain 
can be understood from the fact that the magnetic field 
from the ^^C nucleus splits the mg = 1 and ms — — 1 
sublevels of the NV~ center, making them insensitive, to 
first order, to the strain. A similar immunity to strain is 
observed when an external magnetic field is applied. 
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At high microwave powers, one can also see resonances 
in group B with a separation of As = 36(2) MHz, group 
D with a separation of A^i = 258(2) MHz, and group E 
with a separation of — 295(2) MHz. Side resonance 
groups B, D, and E are symmetrically displaced from 
the center resonance. As seen in Fig. [2(b)] and Table |l) 
the contrast of resonance groups D and E requires higher 
microwave powers to saturate than the resonance groups 
A and C. This implies that the D and E resonances are 
of different origins than the C resonances. 



A. Group C resonances 

By measuring ODMR signals of samples with different 
concentrations of substitutional nitrogen and equal 
concentrations of ^'^C, we can gain information about 
the relevance of substitutional nitrogen and ^^C to the 
signals. ODMR signals taken with the W2 and W5 
diamond samples are shown in Fig. |4] The essential 
differences between these two samples are the concen- 
trations of substitutional nitrogen and NV^ centers. 
The amount of ^'^C in the two samples is the same; it 
is set by the natural abundance of ^^C of 1.1%. Side- 
resonance group C is present in both the signals of the 
low-nitrogen-concentration sample W2 and those of the 
high-nitrogen-concentration sample W5, while the side 
resonance groups B, D, and E are absent in the signals of 
the W2 sample. Since the concentration of ^^C does not 
vary across samples, the resonances in group C are most 
likely due to interactions of the NV~ center with ^'^C. 
Moreover, this matches the previously-observed side 
resonances with separation of 130 MHz attributed to 
hyperfine coupling of the NV~ center and a ^'^C nuclear 
gpjjJSIJel^ The origin of the other side-resonances possibly 
involves coupling of the NV~ center to substitutional 
nitrogen. 

Hyperfine coupling to C nuclear spin 

Besides presenting the experimental data, we also pro- 
vide detailed theoretical calculations which predict the 
positions of all the side resonances. Here, we give a brief 
discussion of the calculations of relevance for the side 
resonance group C, for more details see Sec. |V A| 

As mentioned before, diamond consist of 1.1% ^'^C 
which has nuclear spin 1/2. When a ^^C atom is located 
as a nearest-neighbor to the vacancy, the hyperfine inter- 
action between the NV~ center electronic spin and the 
^■^C nuclear spin is particulary strong. The hyperfine cou- 
pling can be written as "Hhf = S • A • I, where here S is the 
NV~ center electronic spin, I is the ^'^C nucle ar sp in, and 
A is the hyperfine tensor. As detailed in Sec. V A and il- 
lustrated in Fig. |3(a) we find that the energy eigen-states 
shift due to the hyperfine interaction. Using known val- 
ues for the components of the hyperfine tensor, we calcu- 
late that magnetic resonance should occur at frequencies 
(2870 - 60) = 2810 MHz and (2870 -I- 69) = 2939 MHz. 
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FIG. 4: Comparison of ODMR signals of CVD and HPHT di- 
amonds. The nitrogen concentration is different between the 
two samples, while the ^^C concentration is the same. Mea- 
surements were done at zero magnetic field with a laser power 
of 1 mW, on diamond samples W2 and W5. The resonant mi- 
crowave Rabi frequency was 1.8 MHz. The W2-sample signal 
is offset by 0.01 for clarity. 



The separation between the side resonances is calculated 
to be 129 MHz. In an ensemble measurement, where only 
a fraction of the NV~ centers have a nearest-neighbor ^^C 
atom, we predict three magnetic resonances; one strong 
resonance with frequency around 2870 MHz and a pair 
of weaker side-resonances separated by 129 MHz. We 
finally note that the measured positions of the C side 
resonances agree with the calculated values to within a 
few MHz, which supports that the C side resonances are 
related to the hyperfine interaction with ^-^C. 



B. Group B, D, and E resonances 

In the presence of a magnetic field, the nominally de- 
generate \ms = ±1) states split according to the Zee- 
man effect. The splitting is illustrated in Fig. 1(a) 
and is equal to 2'-fBz, where is the projection of 
the external magnetic field on the NV~ axis and 7 = 
1.761 X lO'^^ s~^T~^ is the gyromagnetic ratio for the 
NV~ center. Magnetic resonance between \ms = 0) and 
\ms — ±1) states therefore occurs at frequencies D±^Bz, 
where D « 2870 MHz is the central resonance frequency 
at zero magnetic field. Note that there exist four pos- 
sible orientations (along the four [111] crystallographic 
directions) of the NV^ center in the diamond lattice. 
Furthermore, the magnitude of the projection B^ on the 
NV" axis depends on the orientation of the NV^ center. 
A significant shift of the side resonances relative to the 
unperturbed NV~ ODMR resonances when a magnetic 
field is applied would suggest that the origin of the side 
resonances is interaction of the NV~ center with another 
system which exhibits the Zeeman effect, i.e. a system 
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FIG. 5: (a) ODMR measurements of sample W5, taken at 
a light power of 1 mW, at room temperature, and at a mi- 
crowave power corresponding to a resonant microwave Rabi 
frequency of 3.0 MHz. A magnetic field was applied along 
the [100] direction; the strength corresponding to each mea- 
surement is written under the curve. Note that the magnetic 
field along the NV~ axis is — B cos 54.5°. The values used 
for labeling purposes reflect the total B value. Signals are 
normalized and offset for clarity. This figure shows some rep- 
resentative spectra of the ones recorded, (b) Solid lines show 
magnetic field dependence of the PI center transition energies. 



Points are extracted from the experimental data in Fig. 5(a) 



Circles and squares on lines of higher and lower frequencies, 
respectively, represent frequency differences between E-group 
and B-group side resonances, respectively, and their associ- 
ated A-group resonances. 



which has an electronic spin (since nuclear spins couple 
weakly to magnetic fields compared to electron spins). 



Figure 5(a) shows ODMR measurements of the W5 



sample in which magnetic field was varied from G 
to 56 G along the [100] crystallographic direction. 
In this direction of applied magnetic field, all four 
orientations of the NV^ center have equal projections 
Bz of the magnetic field on the NV~ axis. Note that 
the diamond sample used has a high concentration of 
nitrogen, and side-resonance groups B, D, and E are 
present. As in Fig. 5(a) in the presence of a magnetic 
field, we observe many more side resonances than at 



zero magnetic field. For simplicity, we will focus on the 
side resonances furthest away (for example, E^) and 
those closest (B^) to the central resonance {A^). The 
outermost side resonances belong to the E-group and 
the innermost side resonances belong to the B-group; 
their offsets from the associated center resonance A are 
plotted in Fig. |5(b)| together with their theoretically 
calculated positions in the presence of a magnetic field. 
The other lines in Fig. |5(b)| correspond to shifts of other 
resonances. The details of the calculations can be found 
in Sections |VB[ |V C[ and |VD| In the absence of an 
external magnetic field, the results of these calculations 
are As « 44 MHz, A^, « 254 MHz, and A^ « 298 MHz, 
which agree with the measured values to within a few 
MHz. The model used for calculating these frequency 
separations has two input parameters which are related 
to the hyperfine interaction between the PI center 
electronic spin (S = 1/2) and the PI center nuclear 
spin (/ = 1). The two parameters are the hyperfine 
interaction strength along the PI center axis (Ay) and 
perpendicular to the PI center axis (^j_). In order to 
calculate the magnetic field dependence, we also assume 
that the PI center is oriented along one of the four [111] 
crystallographic directions. These measurements and 
the agreement of the positions of the side resonances 
with the theory derived in Sections |VB[ |V C[ and |VD| 
further support interaction of the NV^ center with the 
PI center as the origin for the side-resonance groups B, 
D, and E. 

Simultaneous spin flips of NV'^ and PI centers 

Apart from correctly predicting the positions of the 
resonances at various magnitudes of magnetic field, the 



theory discussed in Sections V B V C and V D gives us 



an intuition of why the interaction between PI and NV^ 
centers leads to side resonances in the ODMR spectra. 
We provide a short discussion here. PI centers have elec- 
tronic spin 1/2 and nuclear spin 1. These impurities are 
present in HPHT-synthesized diamonds with concentra- 
tion of approximately 100 ppm, in comparison to the 
^^C concentration of 10'' ppm. Due to the low concen- 
tration of PI centers, the average separation between 
PI and NV~ centers is large, leading to a weak mag- 
netic dipole-dipole coupling. This weak coupling slightly 
mixes the energy levels of the two, which allows simulta- 
neous spin flips driven by the applied microwave field^. 
Such simultaneous spin flips result in side resonances 
displaced from the unperturbed NV~ resonance by an 
amount corresponding to energies of transitions within 
the PI center. In the absence of a magnetic field, the 
PI center has three energy levels (each two- fold degen- 
erate) [see Fig. |3(b)| , leading to three transitions within 
the PI center. The transition energies are calculated in 
detail in Sections |VB[ V C and VD and are found to 



be {22, 127, 149} MHz. This gives rise to three groups of 
sidebands symmetrically located around the unperturbed 
NV" resonance A with separations of {44, 254, 298} MHz 
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as illustrated in Fig. |3(b)[ 



V. APPENDIX: THEORY 



IV. CONCLUSION AND OUTLOOK 

Experimentally and theoretically we identify the ori- 
gins of side resonances to the unperturbed NV~ magnetic 
resonance which at zero magnetic field occurs at approx- 
imately 2870 MHz. We attribute symmetrically located 
side-resonance groups separated by around 40 MHz, 260 
MHz, and 300 MHz to simultaneous flips of spins of the 
NV~ center and weakly coupled PI centers in the dia- 
mond. Our results confirm the assignment of the asym- 
metrically displaced side resonances separated by around 
130 MHz to hyperfine interaction of the NV~ center with 
a nearest-neighbor ^'^C nuclear spin. 

The side resonances due to coupling between NV~ and 
PI centers are of importance for magnetic-field sensing 
with NV~ centers in nitrogen rich diamonds. The ex- 
istence of many side resonances complicates the ODMR 
spectrum and can lead to significant broadening of the 
NV~ magnetic resonances when the side resonances over- 
lap with the unperturbed NV~ resonances. Simultane- 
ous spin-flips can be thought of as a source of relaxation 
which depolarize the NV~ center, leading to a reduced 
longitudinal spin- relaxation (Ti) time. The reduction 
of Ti was observed indirectly in Ref. [T3| as a nonlinear 
change in the linewidth of the NV~ magnetic resonance 
as a function of resonant microwave Rabi frequency. 

The side resonances due to coupling of NV~ and PI 
centers could perhaps be used to polarize and detect PI 
centers, which, in contrast to the NV~ centers, cannot be 
directly optically polarized or probed. By driving several 
side resonances sequentially, it should be possible to fully 
polarize and initialize the PI center in a specific quantum 
state. The ability to polarize PI centers is important for 
magnetometry and quantum information processing with 
NV~ centers. Previously, PI centers have been thermally 
polarized by cooling a diamond located in a high mag- 
netic field down to a few kelvirfi^. This polarization leads 
to improved spin relaxation times Ti and T2, which are 
of importance for magnetic field sensing^^^^. 
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A. Hyperfine coupling to ^ C 

Diamond consists mainly of ^^C atoms which have no 
nuclear spin. However, '^'^C, which is present due to the 
1.1% natural abundance, has nuclear spin 7=1/2 and 
can couple to the NV~ center electronic spin S* = 1 
through the magnetic dipolar (hyperfine) interaction. 
This interaction is particularly strong when the -'^'^C atom 
is nearest-neighbor to the vacancy in the NV~ center. 
The nearest-neighbor hyperfine interaction between the 
NV~ center and a ^'^C atom has been studied using elec- 
tron paramagnetic resonance (EPRJPEl and optically de- 
tected magnetic resonancd^. 

Here, we calculate the energy spectrum of the com- 
bined system consisting of the electronic spin of an 
NV~ center and a nearest-neighbor ^■^C nuclear spin [see 
Fig. 3(a)| . In the absence of strain and magnetic fields, 
the NV~ center has C^y symmetry. The ground states 
\ms — 0) and \ms = ±1) are separated by D « 2.87 GHz. 
In this case, the Hamiltonian for the NV~ center can be 
written as Ho = D{Sz)^- The quantization axis for the 
NV~ center, here denoted z, is along one of the four [111] 
crystallographic directions. 

The hyperfine interaction energy can be written as 



(2) 



where A is the hyperfine tensor. For the case of the NV~ 
center and a nearest-neighbor ^^C atom, the hyperfine 
interaction has axial symmetry around the axis, here de- 
noted by z' , which connects the vacancy and the ^^C 
atonP. I.e., one of the [111] directions which is not the 
NV^ axis. The angle between the z and z' directions is 
e w 109.47°. 

We consider the two coordinate systems X and X' with 
axes x,y,z and x',y',z'. We choose the y-axes to coin- 
cide: y = y' . X' can be obtained from X by a rotation 
around the y-axis with angle 9. We have the relation 



' cos 9 0— sin 

1 
, sin 9 cos 




no) y 



(3) 



We can express the hyperfine Hamiltonian in either co- 
ordinate system 



n 



hf ^ 



S • A • I = S' • i' • I', 



(4) 



where S and I arc the electronic and nuclear spins, re- 
spectively, in the coordinate system X, and S' = TZ{9)S 
and I' = TZ{9)I in X' . 

Since the hyperfine interaction has axial symmetry, the 
hyperfine tensor has a simple form in the coordinate sys- 
tem X': 



A' 



'A^ 
v4_L 
All 



(5) 
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The hyperfine couphng parameters have been measured 
in previous works and have values A\\ 199.7(2) MHz 
and Aj_ « 120.3(2) MRz^. 

We can obtain the hyperfine tensor in the coordinate 
system X using the transformation 



7^(6')U'7^(e') 




(6) 



where we have defined the constants -B|| = y4|| cos^ 9 + 
A±sin^6 and B± = {A± - A||) sin 61 cos 0. Using the 
above-mentioned values for j4||, A±, and 6, we calculate 
B|| = 129.1(2) MHz and B±_ = -25.0(1) MHz. 

We now find the eigen-energies of the combined system. 
Notice that the NV~ zero-field splitting D is much larger 
than the hyperfine-coupling parameters and A±. Us- 
ing the secular approximation, where only the terms in 
Hhf which are proportional to Sz are kept, we find the 
total Hamiltonian 

HIT - Ho + - D{Szf + B^SJ, + B\\SJz. (7) 

The Hamiltonian can be written in matrix form using the 
basis {\ms,mi) \ms = — 1,0, 1, rn/ = —1/2, +1/2}. By 
diagonalizing this Hamiltonian, we find the eigen-energies 
0, D ± {Bf, + ^2)1/2/2, each with a de generacy of 2 [see 



Fig. 3(a) . The combined system of an NV~ center and a 
nearest-neighbor ^'^C atom has two magnetic resonances 
with frequencies D±{B'^^+B\Y/'^ /2. Inserting the values 
for B\\ and Sj^, we find magnetic resonance frequencies 
of 2.87 GHz ± (1/2) • 131.5(2) MHz. Notice that with 
the secular approximation, one calculates that the two 
resonances should be symmetric around the unperturbed 
NV~ resonance at £> « 2.87 GHz. 

In the experiment, we observe that the resonance group 
C due to ^■^C is not symmetrically displaced from the cen- 
tral resonance A [see Fig. 2(a) . This can be explained 



by corrections beyond the secular approximation. We 
can diagonalize the full Hamiltonian Htot = + Hm- 
In this case, we find the resonance frequencies to be 
-60.2(1) MHz and +68.7(1) MHz displaced to the low- 
and high- frequency sides of the central resonance at 
around 2870 MHz. The separation between the two res- 
onances is in this case 128.9(2) MHz. Notice that by 
full diagonalization, we calculate that the resonances are 
asymmetrically located with respect to the unperturbed 
NV^ resonance. 



B. Energy level structure of the PI center 

We here calculate the energy- level structure of the PI 
center in diamond. The PI center has electronic spin 
S = 1/2 and nuclear spin 1=1 which are coupled by the 
electron- nuclear hyperfine interaction "Hm = S- A-I where 
A is the hyperfine-interaction tensor. Since the hyperfine 
interaction of a PI center has axial symmetry along one 



of the four possible [111] crystallographic direction^^H^ 
the hyperfine Hamiltonian can be written as 



n 



hi 



■ S^Iz + • {Sxlx + Syly). 



(8) 



The z-direction is here defined as the direction where the 
PI center hyperfine coupling is largest. We use values 
Aj_ = 81 MHz and A\\ = 114 MHz which has been mea- 
sured using EPIft2^ (uncertainties on A± and A|| are not 
given in Ref. [2^. 

First, we construct the ordered orthonormal tensor 
product basis B in the form [rng,™/) from the electron 
spin basis Bs = {[- 1/2) . |l/2)} and the nuclear spin 
basis B/ = {|-1),|0),|1)},' 



B=< -1,-1 



2' 



-1 



\.o). -i.l), i.l 



Using this ordered basis, we find the following block- 
diagonal matrix: 



n = 



The eigenenergies of this system are 
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(9) 

each with a degeneracy of 2. Plugging in numerical values 
for and A±, we calculate the following values for the 
eigenenergies {-92, 35, 57} MHz. Between these states, 
there are three transitions with frequencies {22, 127, 149} 
MHz. 



C. Simultaneous spin flips of NV and PI centers 

From experiment, we determined that resonance 
groups B, D, and E [see Fig. 2(a)| are due to coupling 
of the NV~ center to PI centers. More specifically, 
these resonances are due to simultaneous spin flips of 
the NV~ center and a PI center. The two impurities 
can interact through the magnetic dipole-dipole interac- 
tion. The interaction strength depends on the distance 
between the NV~ and PI center and is expected to be 
weak (for instance, Ref . IT^ gives 1-10 MHz for the inter- 
action strength). This dipole-dipole interaction allows 
for simultaneous spin-fiips at frequencies corresponding 
to the sums (spin-flip in the same direction) and differ- 
ences (spin-flip in opposite directions) between the typi- 
cal transitions of NV~ centers (on the order of GHz) and 
the possible transitions of the PI centers (on the order 
of MHz) described above. 
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FIG. 6: The theoretically-predicted dependence of the 
eigenenergy levels of the PI center on the magnetic field ap- 
plied in the [100] crystallographic direction is illustrated. 



where (j) is the angle between the direction of the mag- 
netic field and the z-axis, which is defined as the direction 
along which the PI center hyperfine interaction is largest, 
i.e. one of the [111] crystallographic directions. Notice 
that both the hyperfine interaction and the magnetic field 
give rise to preferred directions in space. The PI center 
can be oriented along one of the four possible [111] direc- 
tions. For a general magnetic field direction, the energy 
of the PI center depends on its orientation. However, 
if the magnetic field is oriented along one of the four 
possible [100] crystallographic directions, the PI center 
energies are the same for all PI centers, independent of 
their orientation. This is because the angle between any 
[100] direction and any [111] direction is ^ « 54.74°. 



In the absence of strain and magnetic field, we predict 
resonances in the ODMR signal at frequencies of {2721, 
2743, 2848, 2870, 2892, 2997, 3019} MHz [see Fig. [3(b)] . 
In other words, we predict a central resonance at 2870 
MHz, and three sets of resonances symmetric about the 
center. These groups of side resonances should be sep- 
arated by {44, 254, 298} MHz. We observe these reso- 
nances in Fig. 2 (a) I as groups B, D, and E, respectively. 
Furthermore, in Fig. [4j the resonances predicted to origi- 
nate from interaction of the NV~ center with PI centers 
are not observed in the ODMR signal of the CVD sample. 



Diagonalizing the total Hamiltonian leads to the en- 
ergy levels shown in Fig. |6] At zero magnetic field, there 
are three energy levels which are each two-fold degener- 
ate. The zero-field energies are given by Eq. When 
we consider nonzero magnetic fields, the degeneracies 
break to form six energy levels. At high magnetic field, 
the energy levels split into two branches; corresponding 
to the PI center electronic spin either parallel or anti- 
parallel with the magnetic field. 



D. Magnetic-field dependence 

We now calculate the energy levels of the PI center as 
a function of magnetic field. The zero-field energy levels 
were calculated above using the hyperfine Hamiltonian 
■Hhf — S • A • I. The magnetic field can be described using 
the Hamiltonian Hb — 7B • S, where 7 is the gyromag- 
netic ratio. The total Hamiltonian is 

'Htot — T~Lm + T~Lb 

= A\\ ■ SJ, + ■ iS.,h + Syly) + ■■■ (10) 

+ 715] [cos (0)5, + sin (0)5,], 



The possible transition frequencies can be calculated 
from the energy diagram shown in Fig. [6] With 6 energy 
levels, up to 15 transitions are possible. The 15 possible 
transition frequencies are plotted in Fig. |5(b)| for mag- 
netic field values up to 60 G (the range investigated in 
the experiment). Some transitions may be forbidden; for 
instance, when the magnetic field is large we expect only 
three allowed magnetic dipole transitions (corresponding 
to flipping the PI center electronic spin but conserving 
the PI center nuclear spin). 
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